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<j_) ' Dijet production by almost real photons has been studied at HERA with the ZEUS 

detector. Jets have been identified using the cone algorithm. A cut on x® BS , the fraction 
of the photon energy participating in the production of the two jets of highest transverse 
energy, is used to define cross sections sensitive to the parton distributions in the proton 
and in the photon. The dependence of the dijet cross sections on pseudorapidity has 
been measured for x® BS > 0.75 and x® BS < 0.75. The former is sensitive to the gluon 
momentum density in the proton. The latter is sensitive to the gluon in the photon. The 
cross sections are corrected for detector acceptance and compared to leading order QCD 
calculations. 
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1 Introduction 



At the high energies available at HERA, interactions between almost real photons (of virtuality 
Q 2 ~ 0) and protons produce jets of high transverse energy [p], ||, 0, |, ||. The presence 
of a 'hard' energy scale means that perturbative QCD calculations of event properties can 
be confronted with experiment. At leading order (LO) two processes are responsible for jet 
production. The photon may interact directly with a parton in the proton (Fig. |I|a), or it 
may first fluctuate into an hadronic state (Fig. [l|b). In the first case, known as the direct 
contribution, the full energy of the photon participates in the interaction with a parton in the 
proton and the fraction of the photon momentum (x 7 ) participating in the hard process is 
equal to one. The final state of the direct process consists of two jets, the proton remnant 
and the scattered electron. In the second case, known as the resolved contribution, the photon 
acts as a source of partons which then scatter off partons in the proton and the fraction of the 
photon momentum participating in the hard process is less than one. The final state in this 
case includes a photon remnant, continuing in the original photon direction, in addition to two 
jets, the proton remnant and the scattered electron. At higher orders this simple distinction 
between direct and resolved is no longer precisely defined. 

With a cut on jet transverse energy of > 6 GeV, direct photon events probe the parton 
distributions in the proton down to 10 3 , where x p is the fraction of the proton's 

momentum entering into the hard process. This process is directly sensitive to the gluon 
distribution in the proton, and complements indirect extractions in the same x p range |Bj using 
the measurement of F 2 in deep inelastic scattering (DIS). The x p values sampled by the resolved 
photon contribution are typically higher than those of the direct contribution, lying in a region 
where the proton parton distributions are constrained by DIS data. Resolved photon processes 
are directly sensitive to the photon parton distributions, particularly the gluon distribution in 
the photon, down x 7 ~ 0.06. This is not the case in the measurement of F 2 in 77 interactions 
at e + e~ colliders, where highly virtual photons are used to probe almost real photons J7|. In 
collisions between two almost real photons, however, the measurement of jet cross sections has 
recently shown potential to constrain the gluon distribution in the photon [§. [i| . 

In this paper we separate the direct and resolved photon contributions to jet production, and 
present dijet differential ep cross sections which are sensitive to the gluon distributions in the 
proton and photon respectively. In describing the cross sections to be measured, particular 
attention is paid to the way in which direct and resolved processes are defined. Differential ep 
cross sections are presented as a function of jet pseudorapidity for direct and resolved photon 
processes and compared to available LO QCD calculations. 

2 Definition of cross sections 

We measure the cross section for dijet photoproduction, ep — > ey p — > e X, for events in which 
X contains at least two jets of E^ 1 > 6 GeV. In this experiment, photoproduction events are 
defined by demanding that the electron is scattered at small angles and does not emerge from 
the beam pipe. This requirement corresponds approximately to a cut of Q 2 < 4 GeV 2 , giving 
a median Q 2 of ~ 10~ 3 GeV 2 [|]]. The cross section measured is da/dfj, where fj = \{rii +772) is 
the average pseudorapidityQ of the two jets of highest transverse energy, with the requirement 

1 The z axis is defined to lie along the proton direction, and 77 = — ln(tan |) where 9 is the angle between the 
jet and the z axis. 
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that IA77I = \r]i — 772 1 < 0.5. The cross section is measured for 0.2 < y < 0.8, where y = E 1 /E e 
is the fraction of the initial electron energy (E e ) carried by the almost real photon with energy 
Ej. This y interval corresponds to 7^ centre-of-mass energies {W 1P ) in the range 132 GeV 
< W 1P < 265 GeV. 

The possibility of experimentally separating samples of direct and resolved photon events was 
demonstrated in || . However, as the simple definition of resolved and direct photoproduction is 
only unambiguous at leading order, it is important to find a definition which is both calculable 
to all orders and measurable. 

For two-to-two parton scattering in LO QCD, energy and momentum conservation give the 
fraction of the photon energy involved in the hard scatter as 



x 



LO 



Erp. 
partons J - J T 



parton _„parton 



2yE e 



(1) 



where yE e is the initial photon energy and the sum is over the two final state partons. For direct 
photon events, x^° = 1. Since it is not possible to measure partons, we define an observable 
in terms of jets which is analogous to x L °. This observable, called x° BS , is the fraction of the 
photon's momentum participating in the production of the two highest E^ 1 jets. The explicit 
definition is, 



Z^jets-^T e 



7 ~ 2yE e ' 



where now the sum runs over the two jets of highest E^ 1 . In the x® BS distribution thus 
obtained, the LO direct and resolved processes populate different regions, with the direct pro- 
cesses concentrated at high values of x° BS . The peak arising from the direct contribution will 
not necessarily lie exactly at x° BS = 1 due to higher order effects and/or hadronisation, but 
will still correspond to the kinematic region where most or all of the energy of the photon is 
available to probe the proton. The relationship between the definitions of x° BS and x^° is 
dependent upon the hadronisation and showering models used to describe the final state and 
the assumptions made in the calculation of the hard process, as x^ becomes ambiguous even 
at next-to- leading order (NLO). For these reasons, we will present hadronic jet cross sections 
and will not use Monte Carlo models to correct back to parton kinematics. For the purposes of 
this paper, the separation between direct and resolved photoproduction is defined by a cut on 
x° BS . Other model- independent definitions for separating resolved and direct photon processes 



have been suggested, often using a cut on the energy in a cone around the photon direction |10 
We have chosen the definition based upon x^ B because it depends only upon the measurement 
of y, E^ 1 and the 7/'s of the jets, without the need to introduce further variables. 

The I Ar]\ cut ensures that for the fj bins in which the cross sections are measured, both jets lie 
well within the acceptance of the ZEUS detector, without imposing additional cuts at high and 
low if et . Applying this cut has two further benefits, which can be seen as follows. When both 
jets are at equal r] and equal E^ 1 , V§ = 2E^ 1 and so the range of possible choices of scale in 
theoretical calculations is reduced. More importantly, by rewriting the expression for x° BS in 
terms of \Ar]\ and f/, and assuming the jets to have equal transverse energy, one obtains 
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Ir ^ cosh ^. (3) 

A similar expression can be written for the proton, 



xr ^ cosh ^ (4) 

where E p is the incident proton energy. When the jets are at equal pseudorapidities the hy- 
perbolic cosine term takes its minimum value of unity. Thus the minimum available x values 
are probed for a given E^ 1 and fj, and there is a strong correlation between fj and x® BS in the 
direct cross section and between fj and yx® BS in the resolved cross section ]I1J. Applying the 
cut |A?7| < 0.5 brings us close to this situation. 



3 The ZEUS detector and beam conditions 

Details of the ZEUS detector have been described elsewhere ||12| . The primary components 
used in this analysis are the calorimeter and the tracking detectors. The uranium-scintillator 
calorimeter fl3| covers about 99.7% of the total solid angle and is subdivided into electromag- 
netic and hadronic sections with cell size, respectively, of 5 x 20 cm 2 (10 x 20 cm 2 in the rear 
calorimeter i.e. the electron direction), and 20 x 20 cm 2 . The calorimeter has an equal response 
to electrons and hadrons within 3% and the energy resolution achieved in test beams, in terms 
of the energy E in GeV, is a/E = 18 %/y/E for electrons and a/E =35 %/y/E for hadrons. 
The timing resolution of a calorimeter cell is better than at =1.5 /y/E © 0.5 ns. 

The tracking system consists of a vertex detector (VXD) []rjj and a central tracking chamber 
(CTD) |L5| enclosed in a 1.43 T solenoidal magnetic field. The interaction vertex is measured 
with a resolution along (transverse to) the beam direction of 0.4 (0.1) cm. 

To allow a precise measurement of the luminosity via the electron-proton Bremsstrahlung pro- 
cess, electron and photon lead-scintillator calorimeters have been installed inside the HERA 
tunnel, subtending small angles from the interaction vertex |16|]. The small-angle electron 
calorimeter is also used to tag a subsample of photoproduction events in the approximate 
range 10~ 8 GeV 2 < Q 2 < 10~ 2 GeV 2 . 

In 1993 HERA operated with 84 colliding bunches of 820 GeV protons and 26.7 GeV electrons 
with typical beam currents around 10 mA and a luminosity of 0.6 x 10 30 cm" 2 s _1 . From these 
colliding bunches ZEUS collected a total integrated luminosity of approximately 0.55 pb _1 . 
Additional unpaired electron and proton bunches circulated to allow monitoring of background 
from beam-gas interactions. 



3 



4 Data selection and jet finding 



ZEUS data acquisition uses a three level trigger system. During the 1993 data taking period, 
we have selected events which were triggered at the first level on regional or transverse energy 
sums in the uranium calorimeter. Events were also triggered on a coincidence of an electron 
measured in the small-angle electron calorimeter and an energy deposit in the rear calorimeter. 
These events were used to check the efficiency of the calorimeter-only triggers used in the cross 
section calculation. 

At the second level trigger, cuts on calorimeter timing were used to remove events caused 
by interactions between the proton beam and residual gas in the beam pipe in front of the 
calorimeter ||17||. The efficiency of the combination of the first and second level triggers to 



select events in the kinematic region of our cross sections has been determined to be greater 
than 98%. 

At the third level trigger, tighter timing cuts were made to reject events arising from proton 
beam-gas interactions. Cosmic ray events were vetoed using information from the tracking 
chambers and calorimeter. Events with no vertex found by the central tracking chambers, or 
with a vertex found at z < —75 cm were rejected. The measured transverse energy outside 
a 10° cone around the forward beampipe (E^ ne ) was calculated, making use of the vertex 
information from the tracking detectors and associating the energy in a cell with the geometric 
centre of that cell. Events with E^ ne < 12 GeV were rejected. After these triggers, 470,000 
events remained. 

For the final analysis, more stringent cuts using calorimeter timing and tracking information 
are made to further reduce the background from cosmic rays and beam-gas interactions. Two 
additional cuts are made, based upon different measurements of y M: 



1. Events with an electron candidate in the uranium calorimeter are removed if E' e is greater 
than 5 GeV, and if the electron gives a measured y e = 1 — ^-(1 — cos9' e ) < 0.7, where 
E' e and 6' e are the energy and angle of the scattered electron. 

2. A cut is made on the Jacquet-Blondel measurement of y, yjB = J2i(Ei — E zi )/2E e , where 
Egi = EiCos(8i), and Ei and 0, are the energy and polar angle of the calorimeter cell. 
The angle is determined using the measured z- vertex of the event. It is assumed that 
the scattered electron is not seen in the uranium calorimeter and so the sum runs over 
all calorimeter cells. For any remaining events for which the scattered electron did enter 
the uranium calorimeter and either was not identified or gave y e above 0.7, the value 
of yjB under this assumption will be near one. Proton beam-gas events will have low 
values oiyjB- To further reduce contamination from both these sources, we demand that 
0.15 < yjB < 0.7, which is an estimator of the actual y interval of 0.2 < y < 0.8, as 
determined from studies of the energy loss in inactive material in front of the uranium 
calorimeter. 



For the surviving events, jet finding is performed on all cells of the uranium calorimeter using 
a cone algorithm ]nj. The cone radius of the jet finding algorithm is defined such that 
R = J A7^ eH + A<pl ell = 1. The algorithm searches in pseudorapidity-azimuth (r] ce u-(j) ce u) space 
for the cone containing the highest summed Et, removes the cells in this cone and continues 
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the search for the next highest cone. If the summed transverse energy of the calorimeter 
cells within a cone (E^ 1 ) is greater than 5 GeV, and if the Ex weighted pseudorapidity of the 
centre of the cone lies within —1.125 < r] cal < 1.875, then the cells inside this cone are selected 
as a jet. The E^ 1 threshold of 5 GeV is to compensate for energy losses from inactive material 
in front of the calorimeter. This energy loss is corrected for in the determination of the cross 
sections. This rj range is dictated by the available statistics and by the necessity to remain in 
the central region of the calorimeter. Events are selected if they have at least two jets satisfying 
these cuts. 

After these cuts a sample of 12,100 photoproduction events remains. This is reduced to a final 
sample of 4,000 events after a cut on the value of At] reconstructed from the two jets of highest 
transverse energy found in the calorimeter (\Ai] cal \ < 0.5). The backgrounds from cosmic rays 
and beam-gas interactions are 0.6% and 1.0%, respectively, as determined from the number of 
surviving events originating from empty bunch crossings and unpaired proton bunches. The 
contamination from events with a scattered electron in the uranium calorimeter is found to 
be around 2.5% from studies using DIS Monte Carlo events. Studies using the independently 
selected sample for which the scattered electron is seen in the small angle electron calorimeter 
and for which lower energy thresholds in the uranium calorimeter are applied, show that the 
combined efficiency of the first level trigger and the E^^ cut is greater than 90% for events 
with two jets of E^ 1 > 5 GeV. 



5 Resolution of kinematic variables 

In order to estimate how well jets and energies are reconstructed in the ZEUS detector, and 
to study the efficiency of the data selection cuts, we have used the HERWIG 5.7 [19| and 



PYTHIA 5.6 |2(J event generators in conjunction with a detailed simulation of the ZEUS 
detector and triggers. The simulated integrated luminosity used is greater than that of the 
data. As shown in |3J] , these generators are able to provide a reasonable description of the data 
with the minimum pT of a hard scatter set to 2.5 GeV and this value was used in the current 
implementations. The parton distribution sets used were MRSD_ |^TJ for the proton and GRV 
LO for the photon ||22|| . For each generator (HERWIG and PYTHIA), samples of direct and 
resolved photoproduction events were combined according to the generated cross sections. The 
Monte Carlo simulations are only used to correct for detector acceptance and smearing. 

The experimental shifts and resolutions predicted by the Monte Carlo simulations were obtained 
by comparing the true variables with the reconstructed values in the simulated detector. The 
measured value of x® BS is called x^ al and is calculated by inserting E™ 1 , r] cal , and yjs into 
Equation (0). For the jet variables (including x° BS ), the true variables are calculated by 
performing jet finding on the final state particles generated by the HERWIG or PYTHIA 
program before detector simulation. The reconstructed value of y {jjjb) is systematically lower 
than the true value by an average of 0.08, and has a resolution of 0.06. The measurements of 
rj, fj, Arj and x° B exhibit no systematic shift and have resolutions of 0.06 (0.07), 0.05 (0.06), 
0.10 (0.13) and 0.06 (0.07) respectively, for HERWIG (PYTHIA). The reconstructed value of 
E J T et is systematically lower than the true value by 15% (16%) for the HERWIG (PYTHIA) 
simulation with a resolution of 12%. 

The description of the energy response of the uranium calorimeter in the Monte Carlo simulation 
has been checked for events with an electron measured in the small-angle electron calorimeter. 
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The photon energy calculated from E y = E e — E' e was compared with the value reconstructed 
from the calorimeter variable ujb (E° al = yjBE e ). In addition, for jets in the central region of 
the calorimeter, the measurement of the transverse energy available from the ZEUS tracking 
detectors has been used to check the simulation of the calorimeter energy response |J. The 
transverse energy of jets in the forward region of the calorimeter (outside the acceptance of the 
ZEUS central tracking detector) has been compared with the Etj? 1 of these central jets. From 
these investigations we conclude that the Monte Carlo description of the energy response of 
the calorimeter to jets of hadrons is correct to within a possible overestimation of the energy 
of 5%. 

The effects of discrepancies between the Monte Carlo simulations and the data on the measure- 
ment of the cross sections are estimated by varying the selection cuts made on the reconstructed 
kinematic variables, both in the data and the simulated sample and are included in the sys- 
tematic errors on the final cross section values. 



6 Results and discussion 

Fig. 0a shows the distribution of x^ al for the final data sample entering into the cross section 
measurements. There is a clear peak at high values of x^ al which indicates the presence of 
direct type processes. The shape of this uncorrected x° BS distribution and the position of the 
peak are insensitive to the problem of calorimeter energy response because of the presence of 
energy terms in both the numerator and denominator of Equation (0) . The distributions from 
PYTHIA and HERWIG are also shown where the Monte Carlo curves have been normalised to 
fit the direct peak in the data. Although the shape of the direct peak is reasonably reproduced, 
the Monte Carlo simulations fail to describe the rise seen in the data at low x^. This effect was 
also seen in |3J]. In the same Fig., the LO direct contribution from HERWIG is shown separately, 
and indicates that defining the direct process with a cut on x° BS of 0.75 corresponds well to 
the LO definition of direct photon processes as implemented in the simulations. 

Fig. [|b shows the uncorrected transverse energy flow per jet l/NdE^/ddr] around the jet axis 
for events with x^ al > 0.75. For this class of events, both HERWIG and PYTHIA reproduce the 
data distribution well. The same distribution is shown in Fig. [|c for events with x c ® l < 0.75. In 
this case both simulations fail to describe the transverse energy flow in the forward region, as 
was also observed for the jets in the inclusive jet cross sections in || and ||. 

6.1 Direct photon dijet cross section 

The direct photon cross section was evaluated from our data by applying acceptance correction 
factors obtained using the Monte Carlo simulations. Events were selected by cuts on the 
reconstructed kinematic quantities as outlined in Section 4 and a cut of x™ l > 0.75. For these 
cuts, we have evaluated, bin by bin, the acceptance correction for the measurement of the 
cross section da/dfj as defined by the kinematic variables described in Section 2 with a cut of 
x OBSy o.75. The efficiency and purity are evaluated for bins of width 0.25, with centres in the 
range —0.75 < fj < 1.0, using the simulation of the ZEUS detector in conjunction with events 
from the HERWIG and PYTHIA generators. 
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Energy Scale Uncertainty 




(nb) 


Error (nb) 


Error (nb) 


(nb) 


-0.75 


0.63 


0.08 


0.16 


0.08 


-0.50 


1.31 


0.13 


0.10 


0.19 


-0.25 


1.32 


0.14 


0.08 


0.20 


0.00 


1.37 


0.11 


0.06 


0.24 


0.25 


1.29 


0.11 


0.17 


0.15 


0.50 


0.88 


0.10 


0.17 


0.16 


0.75 


0.50 


0.06 


0.27 


0.06 


1.00 


0.23 


0.04 


0.17 


0.05 



Table 1: da/dfj for ep -> eX + 2 (or more) jets, | Arj\ < 0.5, E 3 T et > 6 GeV, 0.2 < y < 0.8, 

Q 2 < 4 GeV 2 , x° BS > 0.75. 



The efficiency is a slowly varying function of fj and is around 50%, falling to 35% in the 
lowest fj bin. The purity in each bin is approximately independent of fj and around 60% and 
depends mostly on migrations across the E^ 1 cut. Taking the ratio of these gives an acceptance 
correction factor which averages around 1.3 and rises to around 1.5 in the lowest fj bin. This 
correction accounts for all detector effects and migrations. We correct back to the final state 
particles and no subtraction of jet pedestal energy {i.e. possible transverse energy around the 
jet direction which is not associated with the hard subprocess) is performed. The differential 
cross section da/dfj for x° BS > 0.75 and the kinematic range indicated in the figure caption, 
is shown in Table [I] and in Fig. |^ where the cross section value is plotted in the centre of the 
bins. 

The systematic uncertainties in the measurement of the differential cross section have been 
estimated by repeating the acceptance correction with both HERWIG and PYTHIA, and by 
varying the cuts made on the reconstructed quantities, and by using different parton distribution 
sets in the simulations. In addition, the bin by bin correction procedure was checked against 
the results of an unfolding method based upon Bayes' theorem f23f , which showed essentially 



the same result. The largest shift in each bin from these variations was found to be similar to 
the statistical error on the data, and is taken to be the total systematic error. 

The systematic uncertainty arising from a possible 5% uncertainty in the mean energies mea- 
sured by the calorimeter is highly correlated between bins, and is therefore excluded from the 
systematic errors and shown separately as a shaded band in Fig. |3|a. Also included in this 
shaded band is the uncertainty in the measurement of the integrated luminosity of 3.3%. The 
shaded band represents the width of the uncertainty around each data point. 

The measured direct dijet cross section is shown in Fig. |^a. The cross section is around 1.3 nb 
at negative fj values, and exhibits a sharp drop near fj = which arises from the cutoff on 
the minimum E^: 1 and the cuts on y. The cross section is compared to several LO QCD 
calculations in which the two final state partons are considered to be jets. The shape of the 
direct cross section differs from that of the LO QCD calculations. However, several effects must 
be considered when comparing data and theory: 

• The gluon distribution of the proton in this kinematic region is not well known. In 
Fig. ||a the data are compared with LO cross section curves [|ll[] calculated using the 
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GRV LO, CTEQ2M @, MRSA ^§ and MRSDg parton distribution sets for the proton. 



The MRSA and CTEQ2M parton distribution sets are global fits to data which include 
HERA measurements of the structure function F 2 ||26||. The contribution from the tail 



of the LO resolved cross section with xtj > 0.75 is included in the calculations using the 
GS2 |27[] photon parton distribution set. This contribution to the LO cross section is small 
(da/dfj ~ 0.1 nb). This is also the case using other available photon parton distribution 



sets except for the LAC3 set pgj] , where due to the high gluon density at high x^ the 
cross section becomes large enough to describe the whole high x° BS cross section without 
any direct component. This parton distribution set is disfavoured by results presented in 
this paper (see below) and by other measurements 0. 

QCD calculations of this cross section are only available at LO. Some estimate of the size 
of higher order corrections can be made by comparing the LO calculation using the GRV 
LO parton distribution set with the other LO calculations using NLO parton distribution 
sets, which in general are around 20% lower[|. In addition, using 2E^ 1 for the hard scale 
instead of E^/2 lowers the cross section by up to 20%. When both jets are at equal r] 
and equal E^, = 2Ej? and so the most common choices of scale are covered by the 
range E 3 T et /2 to 2E J T et . 

Due to the fact that we are probing low x partons in the proton, the standard approx- 
imation that the incoming partons are collinear and on-shell may be invalid. Analytic 
estimates of the effect of allowing the incoming partons to develop non-zero transverse 



momentum (fc^) have been made in |^S] using a so-called 'fc^ factorization' prescrip- 
tion. Calculations of the cross section using this technique are so far available only for 
the photon-gluon fusion contribution to the direct photon cross section. The typical k^ 
developed is of the order of a GeV. In Fig. we show again the data of Fig. |3]a, and 
the standard LO QCD calculation using the GRV LO proton parton distributions. Also 
shown are the curves for the same parton distribution set for the x^° = 1 contribution 



alone, and for just the gluon-induced part of the x^ = 1 cross section. (From this it can 
be seen that according to LO QCD, most of the direct photon cross section is attributable 
to the photon-gluon fusion diagram.) This latter curve may then be compared to the final 
curve, which shows the l k^ result for the same cross section. Non-zero parton kx in this 
prescription lowers the cross section by as much as 30%, and could bring a full calculation 
into better agreement with the data. 

Non-perturbative 'hadronisation' effects can be expected to be significant for jets. In 
Fig. |3|c we show the data of Fig. |3|a compared to Monte Carlo estimates of the cross 
section calculated from partons generated with HERWIG using the GRV LO proton 
parton distribution set and the LAC1 [^8| photon parton distribution set. Also shown is 
the cross section obtained by performing jet finding on the simulated final state particles. 
This histogram agrees with the shape of the data better than the analytic LO QCD 
calculations of Fig. |]a. Over most of the fj region these histograms lie within 20% of each 
other. 



Considering the uncertainties in LO calculations arising from the choice of the hard scale (taken 
to be E T e /2 in these calculations) and the fact that hadronisation effects are not included 

2 The curve calculated with the GRV NLO parton distribution set (not shown) is also below the GRV LO 
curve by around 20%. 
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in the theoretical curves, the description of the measured points by the theoretical curves is 
reasonable. The effects of kr, hadronisation, and of using different parton distributions are of 
comparable magnitude. That the direct jet cross section is sensitive to the parton densities of 
the proton can be seen by comparing different Monte Carlo estimates of this curve generated 
using different proton parton distribution sets. In Fig. |3]d, we show the data compared to 
hadronic jet cross sections estimated with HERWIG using the GRV LO and MRSD proton 
parton distribution sets. The LAC1 photon parton distribution set was used in both cases 
for the LO resolved x° BS > 0.75 contribution. The separation between parton distribution 
sets remains after hadronisation, and the MRSD parton distribution gives a consistently lower 
cross section for negative values of fj. 

In summary, the LO QCD predictions for the direct cross section are consistent with the data 
at the level of 30%. The theoretical cross section is sensitive to the choice of proton parton 
distribution function. Complete NLO calculations will reduce the ambiguity in comparisons to 
the data. 



6.2 Resolved photon dijet cross section 

The acceptance correction for the resolved photon process as defined in Section 2 with a cut 
of x° BS < 0.75 has been evaluated in the same way as in the direct photon cross section 
measurement. The efficiency and purity are now evaluated for bins with centres in the range 
0.0 < fj < 1.5. The efficiency and purity are both approximately flat and around 40%, giving 
an acceptance correction factor of around unity across the whole range of fj. As for the direct 
measurement, the purity depends mostly on migrations across the E^ 1 cut. No subtraction 
of jet pedestal energy is carried out. However, variations in the kinematic selection of the 
analysis now give systematic variations in the cross section which are larger than the statistical 
errors, and thus the systematic error bars are correspondingly larger. As was seen in Fig. |2|c, 
both Monte Carlo simulations fail to describe the forward region of the transverse energy flow 
around the jet axis, and a reduction of these systematic uncertainties will require improvements 
in the simulations used. In Fig. ^ the correlated uncertainty arising from the description of the 
calorimeter energy response in the Monte Carlo is again shown separately by the shaded band. 
The band also includes the 3.3% uncertainty in the measurement of the integrated luminosity. 
The shaded band represents the width of the uncertainty around each data point. 

Table [2] and Fig. |3] shows the measured cross section da/dfj for x° BS < 0.75 and the kinematic 
range indicated in the figure caption. The cross section is around 2 nb for central values of fj 
and rises to 4 nb at fj = 1.5. The LO cross sections, shown for comparison, are calculated using 
different photon parton distribution sets and the MRSA set for the proton. The theoretical 
sensitivity to the parton distributions in the proton is small (not shown), with the variations 
between curves calculated using different parton distribution sets being much less than the 
estimated errors on the measured cross section. As x p is higher in these events than in the 
direct, the effect of kx from partons in the proton on this calculation is expected to be small. 



The DG GRV LO and GS2 parton distribution sets reproduce the shape of the cross section 



well and can be brought into agreement with the data by applying a multiplicative factor of 1.5 



to 2. (Note that the available NLO calculations [[JT[] for inclusive jet photoproduction, which is 
dominated by the resolved process, differ from LO calculations by a factor of up to two.) The 
LO calculations using the LAC1 and LAC3 parton distribution sets cannot be brought into 
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Table 2: da/dfj for ep -> eX + 2 (or more) jets, | Ar]\ < 0.5, E$r > 6 GeV, 0.2 < y < 0.8, 

Q 2 < 4 GeV 2 , x° BS < 0.75. 

agreement with the data by a constant normalisation factor. However, since the jet pedestal 
energies are not well described by the Monte Carlo for the resolved process, it is difficult to 
estimate what the effects of hadronisation and parton showering might be on the shapes of the 
curves, and so no parton distribution set can be completely excluded at this stage, although 
LAC3 is disfavoured. 



7 Conclusions 

Differential dijet cross sections have been measured in photoproduction with the ZEUS detector 
at HERA. The cross sections measured are da/dfj, \Ar]\ < 0.5, E 3 T et > 6 GeV and 0.2 < 
y < 0.8, in the regions of x° BS > 0.75 (direct photoproduction) and x° BS < 0.75 (resolved 
photoproduction). The measured cross sections have been defined in such a way that they are 
calculable to higher orders in QCD. We have corrected back to the final state particles and no 
subtraction of jet pedestal energies has been carried out. Our results are compared with the 
expectations of LO QCD. The direct cross section, which is sensitive to the gluon content of 
the proton, is consistent with LO QCD calculations to within 30%. Hadronisation, incoming 
parton kx and the choice of scale also influence this comparison. The shape of the resolved 
cross section, which is sensitive to the gluon content of the photon, is described by LO QCD 
calculations using the DG, GRV LO and GS2 photon parton distributions. However, most LO 
QCD calculations lie below the data by a factor of 1.5 to 2. Comparison with NLO calculations 
for both the direct and resolved photon cross sections would be extremely valuable and should 
allow stronger conclusions to be drawn. 
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Figure 1: Examples of leading order diagrams for a) direct and b) resolved photoproduction. 
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Figure 2: a) The x^ al distribution. The solid circles are uncorrected ZEUS data. The solid 
(dashed) line represents the distribution from the PYTHIA (HERWIG) simulation. The LO 
direct contribution to the HERWIG distribution is shown by the shaded histogram. The Monte 
Carlo curves have been normalised to fit the direct peak in the data, b) and c) show the 
uncorrected transverse energy flow 1 / TV dE? / 'dSrj around the jet axis, for cells within one radian 
in of the jet axis, for b) direct and c) resolved events. The solid (dashed) line represents the 
distribution from PYTHIA (HERWIG). 
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c) d) 
Figure 3: da/di] for ep — > eX + 2 (or more) jets, \At)\ < 0.5, E 3 ^ > 6 GeV, 0.2 < y < 
0.8, Q 2 < 4 GeV 2 , x® BS > 0.75. The solid circles are corrected ZEUS data. The inner error bars 
indicate the statistical errors, the outer error bars show the systematic uncertainty (excluding 
the correlated uncertainty) added in quadrature. The shaded band shows the correlated uncer- 
tainty from measurement of energy in the calorimeter and the integrated luminosity. In a) the 
data are compared to LO QCD calculations using several parton distribution sets for the proton 
and the GS2 set for the photon. In b) the data are compared to the LO QCD calculation of 
a) using the GRV LO parton distribution set for the proton, the same calculation but omitting 
the resolved contribution with x^°> 0.75, including only gluon induced direct photon processes 
and the { kx factorized' curve. In c) the data are compared to HERWIG Monte Carlo estimates 
of the cross section using partons and final state jets. For these HERWIG histograms we have 
used the the GRV LO (LAC1) proton (photon) parton distribution set. In d) the data are 
compared to HERWIG jet cross sections using the GRV (LAC1) and MRSD (LAC1) proton 
(photon) parton distribution sets. 15 



ZEUS 1993 




Figure 4: da/dfj for ep — > eX + 2 (or more) jets, |A^| < 0.5, > 6 GeV, 0.2 < y < 
0.8, Q 2 < 4 GeV 2 , x° BS < 0.75. The solid circles are corrected ZEUS data. The inner error bars 
indicate the statistical errors, the outer error bars show the systematic uncertainty (excluding 
the correlated uncertainty) added in quadrature. The shaded band shows the correlated un- 
certainty from measurement of energy in the calorimeter and the integrated luminosity. Also 
shown are LO QCD calculations. The parton distribution sets used for the photon are LAC3, 
GS2, GRV, LAC1 and DG. The proton parton distribution set used is the MRSA set. 
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